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ABSTRACT 

The centers of ellipticals and bulges are formed dissipationally, via gas inflows over short 
timescales - the "starburst" mode of star formation. Recent work has shown that the surface 
brightness profiles, kinematics, and stellar populations of spheroids can be used to separate 
the dissipational component from the dissipationless "envelope" made up of stars formed over 
more extended histories in separate objects, and violently assembled in mergers. Given high- 
resolution, detailed observations of these "burst relic" components of ellipticals (specifically 
their stellar mass surface density profiles), together with the simple assumptions that some 
form of the Kennicutt-Schmidt law holds and that the burst was indeed a dissipational, gas- 
rich event, we show that it is possible to invert the observed profiles and obtain the time and 
space-dependent star formation history of each burst. We perform this exercise using a large 
sample of well-studied spheroids, which have also been used to calibrate estimates of the 
"burst relic" populations. We show that the implied bursts scale in magnitude, mass, and peak 
star formation rate with galaxy mass in a simple manner, and provide fits for these correlations. 
The typical burst mass Mburst is ^ 10% of the total spheroid mass; the characteristic starburst 
timescale implied is a nearly galaxy-mass independent fburst ^ 10^ yr; the peak SFR of the 
burst is Mburst /fburst; and bursts decay subsequently in power-law fashion as cx t~^'^. 
As a function of time, we obtain the spatial size of the starburst; burst sizes at peak activity 
scale with burst mass in a manner similar to the observed spheroid size-mass relation, but 
are smaller than the full galaxy size by a factor ~ 10; the size grows in time as the central, 
most dense regions are more quickly depleted by star formation as /?burst « f"-^. Combined 
with observational measurements of the nuclear stellar population ages of these systems - 
i.e. the distribution of times when these bursts occurred - it is possible to re-construct the 
dissipational burst contribution to the distribution of SFRs and IR luminosity functions and 
luminosity density of the Universe. We do so, and show that these burst luminosity functions 
agree well with the observed IR LFs at the brightest luminosities, at redshifts z ~ — 2. At 
low luminosities, however, bursts are always unimportant; the transition luminosity between 
these regimes increases with redshift from the ULIRG threshold at z ^ to Hyper-LIRG 
thresholds at z 2. At the highest redshifts z > 2, we can set strict upper limits on starburst 
magnitudes, based on the maximum stellar mass remaining at high densities at z = 0, and 
find tension between these and estimated number counts of sub-millimeter galaxies, implying 
that some change in bolometric corrections, the number counts themselves, or the stellar IMF 
may be necessary. At all redshifts, bursts are a small fraction of the total SFR or luminosity 
density, approximately 5 — 10%, in good agreement with estimates of the contribution of 
merger-induced star formation. 

Key words: galaxies; formation — galaxies: evolution — galaxies: active — star formation: 
general — cosmology: theory 



1 INTRODUCTION 
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Understanding the global star-formation history of the Universe re- 
mains an important unresolved goal in cosmology. Of particular in- 
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terest is the role played by mergers in driving star formation and/or 
the infrared luminosities of massive systems. A wide range of ob- 
servations support the view that violent, dissipational events (e.g. 
gas-rich mergers) are important to galaxy evolution, and in partic- 
ular that the central, dense portions of galaxy bulges and spheroids 
must be formed in such events; but less clear is their contribution 
to the global star formation process. 

In the local Universe, the population of star-forming galax- 
ies appears to transition from "quiescent" (undisturbed) disks - 
which dominate the total star formation rate/IR luminosity density 
- at the luminous infrared galaxy (LIRG) threshold lO' ' Lq (M, ~ 
10 — 20 Mq yr~ ' ) to systems that are clearly merging and violently 
disturbed at a few times this luminosity. The most intense starbursts 
at z = 0, ultraluminous infrared galaxies (U LIRGs;Lir > 10''^L,^) , 
are in variably associated with mergers (e.g. [Joseph & Wright! 19851 ; 
I Sand ers & Mirabel 1996), with dense gas in their centers provid- 
ing material to feed black hole growth and to boost the concentra- 
tion and central ph ase space density of merging spirals to matc h 
those of ellipticals l lHemauistet"Z[| 19931 ; iRobertson et"ai]|2006h . 
Various studies have shown that the mass involved in these star- 
burst events is critical to explain the relations between spirals, 
mergers, and ellipticals, and has a dramatic impact on the proper- 
ties of me rger remnan ts (e.j 



.Lake & 

19941; [ShT er & Fisch eill 1998; I James et al.lll999l;lGenzel et al.ll200ll; 



Tacconi et al. 2002; Dasvra et al."2006l. 120071 ; iRothberg & Josephl 
2004. 2006a; Hopkins et al. 2009b,9). 

At high redshifts, bright systems dominate more and 
more of the IR luminosity function (e.g. Le Floc'h et al.l |2005|; 
Perez-Gonzalez et al.l |2005| ; ICaputi et all |2007| ; iMagnelH et al.l 
20091) . Merger rates increase rapidl y (by a factor ~ 10 from z = 



— 2; see e.g. iHopkins et al. I l2009ll , and references therein), lead- 
ing to speculation that the merger rate evolution may in fact drive 
the observed evolution in the cosmic SFR density, which also rises 
sharply in this interval (e.g. lHopkins"& Beacom ,200^, and ref- 
erences therein). However, many LIRGs at z ~ 1, and possibly 
ULIRGs at z ~ 2, appear to be "normal" galaxies, without dra- 
matic morphological disturbances associated with the local star- 
burst population or large apparent AGN contributions ( Yan et al.' 
2007 ; Sajinaetal. 2007; Dey et al. 2008; Melbourne et al. 2008; 
Dasvra et alj2008l) . At the same time, even more luminous systems 



appear, including large numbers of Hyper-LIRG ( HyLIRG; Lir > 
lO'^Lp) and bright sub-millimeter galaxies (e.g. IChapman et al.l 
I2OO5I ; lYounger et alllioOTl ; ICasev et alj|2009l) . These systems ex- 
hibit many of the traits more commonly associated with merger- 
driven starbursts, including morphological disturbances, and may 
be linked to the emergence of massive, quenched (non star- 
forming), compact ellipticals at times as early as z ~ 2 — 4 



("Papovich et al. 2005; Younger et al. 


'2008; Tacconi et al." 200^ 


fSchinnerer et al. 2008; Chapman et al 


2009; Tacconi et al. 200^ 


In a series of papers, Hopkins et alj ( 2008d. l2009blle() (here- 



after H08c, H09b, H09e, respectively), the authors combined simu- 
lation libraries of galaxy mergers with observations of nearby ellip- 
ticals to develop a methodology to empirically separate spheroids 
into their two dominant physical components. First, a dissipation- 
less component - i.e. an "envelope," formed from the violent re- 
laxation/scattering of stars already present in merging stellar disks 
that contribute to the remnant. Because disks are extended, with 
low phase-space density (and collisionless processes cannot raise 
this phase-space density), these stars will necessarily dominate the 
profile at large radii, hence the envelope, with a low central density. 

Second, a dissipational or "burst" component - i.e. a dense 
stellar relic formed from gas which lost its angular momentum and 



fell into the nucleus of the remnant, turning into stars in a compact 
central starburst like those in e.g. local ULIRGs (although these 
represent the most extreme cases). Because gas radiates, it can col- 
lapse to very high densities, and stars formed in a starburst will 
dominate the profile within radii ~ 0.5 — 1 kpc, accounting for the 
high central densities of ellipticals. The gas will reflect that brought 
in from merging disks, but could in principle also be augmented 
by additional cooling o r stellar mass loss in the elliptical (see e.g. 
ICiotti & Ostrikeill2007h . 

In subsequent mergers, these two stellar components will act 
dissipationlessly, but the segregation between the two is suffi- 
cient that they remain distinct even after multiple, major "dry" re- 
mergers: i.e. one can still, in principle, distinguish the dense central 
stellar component that is the remnant of the combined dissipational 
starburst(s) from the less dense outer envelope that is the remnant 
from low-density disk stars. 

In H08c, the methodology for empirically separating 
these components in observed systems is presented, and 
tested this on samples of nearby merger remnants from 
iRothberg & JoseplJ ( l2004h . Comparison with other, indepen- 
dent constraints such as stellar pop ulation synthe sis models 
jTitus et alJl997l ; lReichardt et al .'2001 ;'Michard' 2006) and galaxy 
abundance profiles 1 Foster et al.. .2009. ; .McDermid et al.. .20061 ; 
ISanchez-BMzquez et alj2007h . as well as direct comparison of sim- 
ulations with observed surface brightness profiles, galaxy shapes, 
and kinematics are used to demonstrate that the approach can reli- 
ably extract the dissipational component of the galaxy (see § 12. It . 

If this general scenario is correct, then the empirically identi- 
fied burst relic components in local spheroids represent a novel and 
powerful new constraint on the history and nature of dissipational 
starbursts. In this paper, we present and develop these constraints 
- and in particular, we show that they are not merely integral con- 
straints on the masses and sizes of bursts. We show that the unique 
nature of dissipational star formation, together with the existence 
of some Kennicutt-Schmidt type relation between gas surface den- 
sities and star formation rates, means that the relic profiles can be 
inverted to obtain the full time and radius-dependent star forma- 
tion history of each galaxy starburst (§[2j. We present the resulting, 
derived burst star formation histories for samples of hundreds of lo- 
cal, well-observed galaxy spheroids. We show that such starbursts 
follow broadly similar time-dependent behavior, with characteris- 
tic starburst star formation rates, durations, rise and decay rates, 
and spatial sizes that scale with galaxy mass and other properties 
according to simple scaling laws across ~ 5 decades in starburst 
and spheroid mass (§ 13. 2t . Combining these inferred star formation 
histories (and resulting burst lightcurves) with empirical determina- 
tions of the ages of each starburst, we can reconstruct the starburst 
history of the Universe, including e.g. the luminosity functions of 
such dissipational bursts at all redshifts, and their contribution to 
the global star formation rate density (§ |3.3l l. We discuss the uncer- 
tainties in this approach, compare with the previous, independent 
attempts to constrain these quantities via other observational meth- 
ods, and summarize our results in §|4] 



Through out, we adopt a Q,m ~ 0.3, ~ 0.7, h = 0.7 cos- 
mology and a IChabried W&t) stellar IMF, but these choices do not 
affect our conclusions. 
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2 THE OBSERVATIONS AND METHODOLOGY 
2.1 The Observations 

As discussed in § [T] splieroid mass profiles can be decomposed 
into a central relic starburst and an outer stellar envelope. The two 
components are distinct physically in the sense that the amount and 
distribution of the two are determined by dissipational and dissipa- 
tionless dynamics, respectively. In H08c, H09b,e, the authors com- 

S iile la rge sa mples of elli p ticals from the studies of lKormendv et al] 
200^ and iLauer et al] ( l2007h . and present decompositions for 
these galaxi es. We adopt these resul ts for the present study. 

Brieflv. lKormendv et al.l l [2009h present a V-band Virgo ellipti- 
cal survey, based on the co mplete sample of Virg o galaxies down to 
extremelv faint systems inbinggeli et al.l (l985l) ( the same sample 

alJI 



tion properties - i.e. independent constraints from stellar popu- 



studied in ICote et alj|200^ ; 



Ferrarese et 



2006h . lKormendvetal.l 



eluding iBender et al.l 


1988|; 


Caon et al.l Il990l, 1 19941; iDavis et al.l 


19851; Kormendv et al. 


20051; 


Laueilll985f; Lauer et alj|l995l 20051; 


Liu et alJl2005l;|Peletier et al. 


|l990|) and new photometry from the 



McDonald Observatory, the HST archive, and the SDSS for each of 
their objects which, after careful conversion to a single photometric 
standard, enables accurate surface brightness measurements over 
a wide dynamic range (with an estimated zero-point accuracy of 
±0.04 K mag arcsec"^). Typically, the galaxies in this sample have 
profiles spanning ~ 12 — 15 magnitudes in surface brightness, cor- 
responding to a range of nearly four orders of magnitude in physical 
radii from ~ lOpc to ~ 100 kpc, permitting the best simultaneous 
constraints on the shapes of both the outer and inner profiles of 
any of the objects we study. Unfortunately, since this is restricted 
to Virgo ellipticals, the number of galaxies is limited, especially at 
the intermediate and high end of the mass function. 

We therefore include surface brigh tness profile s from 



Lauer et al 



Lauer et al 



2007 ), further supplemented by iBender et al. Il ll988l) . 
2007n compile V-band measurements of a large num- 



ber of nearby systems for which H ST imagine of ga lactic nu- 
clei is avail able. These include the iLauer et all 1 I2OO5I) WFPC2 
data-set, the iLaine et alj ( l2003ij) WFPC2 BCG sample (in which 
the o bjects are specific ally selected as b rightest clus t er gal axies 
from IPostrnan & Lauerl Jl995i) ). and the iLauer et al] ( 1 19951) and 
iFaber et al. I h997h WFPCl compilations. Details of the treatment 
of the profiles an d conversion to a single standard are given in 
iLauer et alj l l2007h . HST images are combined with ground-based 
measurements (see references above) to construct profiles that typ- 
ically span physical radii from ~ lOpc to ~ 10 — 20 kpc. The sam- 
ple includes ellipticals over a wide range of luminosities, down to 
Mb ~ — 15, but is dominated by intermediate and giant ellipticals 
and SO galaxies, with typical magnitudes Mb J; —18. 

H08c and H09b,e apply various tests to demonstrate that the 
profile decompositions are robust. Systematic comparison with 
large suites of hydrodynamic simulations as calibrators implies 
that the derived burst masses and radial profiles are accurate to 
within a factor of ~ 2. Indepe ndent, purely empirical analyses of 
thes e profiles (K ormendv et al. 200 9i) and those of similar objects 
(e.g. lFerrarese et aUi2006. ; ,Balcells et alj|2007ah yield similar con- 
clusions. For most of the objects in our data-set, the profiles in- 
clude e.g. ellipticity, an /a, and g — z col ors as a f unction of ra- 
diu s, as well as kinematic information. In iKormend v et al. (200§) 
and lHopkins et aL I( l2009hll3) . the authors show that these secondary 
properties exhibit transitions that mirror the fitted decompositions 
(e.g. transitions to diskier, more rotationally supported, younger in- 
ner components). 

Likewise, comparison with detailed resolved stellar popula- 



- yields good agreement where available iMcDermid et al. 2006 


Sanchez-Blazauez et alJl2007;lReda et al.ll2007l;IFoster et al.ll200S 


Schweizei I996I; Titus et alj 19971; 


Schweizer & Seitzed Il998 


20071; iReichardt et al.ll200ll; iMichard 


2006|). Moreover, in voung 



merger remnants, where stellar populations and colors more clearly 
indic ate the post- starburst componen t, the methodology wo rks 
well jRothberg & Josephll2004l2006al lbl; lHoDkins et al.|[2008ch . In 
H09b,e, these results are also checked against profile data used in 
IBender et al] ( Il988lll99lll99l ll9941. The latter are more limited 
in dynamic range, but allow for construction of multi-color profiles 
in e.g. V , R, and / bands, to test whether they are sensitive to the 
band used and to construct e.g. stellar mass and M*/L profiles. 

We illustrate our methodology in Figure [T] First, we show a 
typical stellar mass surface density profile of an ~ L, elliptical 
galaxy, from the sample of iLauer et alj j2007h . We plot the de- 
composition into burst and extended envelope (violently relaxed) 
components, together with some of their integral properties (stellar 
mass fractions, effective radii, and fitted profile shapes in the form 
of the best-fit Sersic index). We compare the observed system to 
the stellar mass density profile of one of the hydrodynamic simula- 
tions presented in H09b. The simulation is determined in that paper 
to be similar in total mass, profile shape, and kinematic properties 
to the observed system (see Table 1 therein). In the simulation, the 
true, physical starburst component from gas that loses angular mo- 
mentum and falls to the center at coalescence is known, as is the 
non-starburst component from violently relaxed stars present in the 
disks before the final merger. We show the separate contributions 
of the two components to the stellar mass density. The agreement 
is good, lending confidence to our empirical de-composition pro- 
cedure, and allowing us to consider the simulations as a reasonable 
calibration sample for our approach. 



2.2 Basic Assumptions 

Taking the observationally estimated burst component, this defines 
a relic surface density Ebursi(^). Of course, in general, the star for- 
mation history that yields a given stellar surface density is non- 
unique. But, if we make three simple, physically and observation- 
ally motivated assumptions, Eburst(^) can be inverted into a space 
and time-dependent burst star formation history. 
These assumptions are: 

• (1) That some form of the lKennicutll i ll 9981) law holds for these 
objects, relating their star formation rate surface density E, to their 
gas surface density Egas. 

• (2) That the relic burst mass is dominated by of order one 
most massive event, i.e. is not constituted from the later assembly 
of many independent, well-separated bursts nor by an extended se- 
ries of independent, well-separated bursts in the same galaxy. There 
may be many small events, but most of the mass should come from 
a single, dominant event. 

• (3) That this event, at its inception, involved the central regions 
being (however briefly) gas-dominated. In other words, the burst 
began from a gas-dominated central density and formed largely in 
situ, rather than from an extended, low-gas-density trickle in which 
the gas density was never nearly as large as the total/final stellar 
mass density. 

In detail, these assumptions are not exactly correct (we discuss 
this further in §|4ll, but most observational constraints indicate that 
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Figure 1. Illustration of the methodology used in this paper. Top Left: Observed stellar mass profile of a typical ~ L* elliptical from the lLauer et alj fcOOTh 
sample (dark blue solid). We show the empirical, fitted decomposition into inner "starburst relic" (dark blue dashed) and outer violently relaxed (dark blue 
dotted) components, from H09b,e. We compare a simulation which yields a remnant with similar profile, kinematics, and stellar populations (total mass profile 
in black). For the simulation, we know the true, dissipational merger-induced starburst component (light blue, with labeled mass fraction), and dissipationless 
violently relaxed component from stars formed before the final merger (red; with labeled Sersic index). The fitting and simulation-matching procedures yield 
similar burst decompositions. Top Right: Implied gas surface density profile of the burst {ignoring the violently relaxed component, formed over much more 
extended periods) as a f unction of time. A t time f = (defined as the time where the implied SFR peaks), we assume the burst relic was in place mostly as 
gas. Given the observed lKennicutj i 19981) relation, this yields a SFR surface density at each location. The gas density then declines as gas is turned to stars. 
Integrating forward, we obtain the gas (and formed stellar) mass distribution at each later time (labeled). Bottom Left: Corresponding SFR surface density 
profile. We also label the half-SFR radius at each time. Bottom Right: Implied total SFR (integrating the SFR density profile over radius), at each time (black 
dashed). The "two-sided" estimate assumes a light curve with a symmetric rise/fall. We compare the true simulation SFR versus time in the best-fit simulation. 
The two agree well - despite the simple nature of the model, it appears to be reliable in simulations with more detailed physics, and allows us to recover the 
star formation history of a given spheroid starburst from its relic mass profile. 



deviations from them are at the factor of ~ a couple level, compa- 
rable to or less than e.g. most of the systematic uncertainties in the 
observations to which we will compare. 

Reg arding assumption (1); observations indicate that at both 
low (e.g. lKenniciit3l 19981 : iBigiel et al.ll2008l : iBothwell et al."2009'. 
and references therein) and high (z ~ 2 — 4, e.g.|Bouche et al. 2007) 
redshifts, a Kennicutt-Schmidt star formation law applies over a 
large dynamic range in stellar mass and star formation rates (from 
~ 1 — 3000 Mq yr^ ' ). These observations have also shown that the 
same relation pertains to a diverse array of objects, from dwarf 
galaxies, through local gas-poor disks, to near "quenched" red sys- 
tems, to clumpy, turbulent high redshift disks and irregular sys- 
tems, to merger-induced starbursts and ULIRGs, and high-redshift 
sub-millimeter galaxies, among others. The extent of the uncer- 
tainty appears to be in the precise normalization/slope of the rela- 
tion (with the references above favoring power-law slopes varying 



from 1 .4 — 1 .7). We take this uncertainty into account below, and 
show that it has a relatively minor effect on our conclusions. The 
lKennicutll ( IT998l) relation may break down at very low surface den- 
sities, where there is some star formation threshold, but our analysis 
here is specifically for very high surface-density systems. 

Regarding assumption (2); this is well-motivated by a com- 
bination of theoretical models and observations of galaxy stellar 
populations and structural properties. In most models, nuclear burst 
masses are dominated by the largest, most rece nt event (see e.g. 
Khoc hfar & Burkert 20 03; C: r oton et al. 2006; de Lucia & Blaizoll 
20071 ; ISomerville et"ar]|2008l ; iHopkins et al J [20091 . Mergers and 
other violent processes that can drive large quantities of gas into 
galaxy centers are not so common as to yield many bursts of equal 
strength, and even if systems had earlier bursts at much higher red- 
shifts, they would have grown by a large factor in mass since those 
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times, such that the most recent event would dominate (se e refer- 
ences above and lWeinzirl et alj2009ljHoDkins et al]|2009lj) . 

Observationally, evidence from nuclear kinematics and pro- 
file shapes supports the view that most of the mass in spheroid 
starburst components, especially in the cusp ellipticals that dom- 
inate (> 90%) the total mass density of spheroids and bulges, was 



cant modification since that time (e.a. Faber et al.l 19971; Kormendv 


19991; 
200 ll; 


Guillen et all 2000l: Rest et all 12001; iRavindranath et al. 


Laine et alj l2003d; iLauer et alj |2005|; Isimien & Prusniel 


1997allbllc; Emsellem et alj 2004, |2007|; iMcDermid et al. 


2006; 


CaoDellari et alj |2007|; iLauer et alJ |2007|; iFerrarese et al. 


2006; 


Cote et al.1 20071; Hookins et alj 2008al; iKormendv et al.1 l2009b. 



Even in core ellipticals, thought to have been affected by subse- 
quent re-mergers, the kinematics require relatively little mixing of 
major populations, with of order a couple similar-mass merged sys- 
tems (i.e. introducing a factor ~ 2 uncertainty) ~ more minor merg- 
ers may be dynamically important, but will (by definition) con- 
tribute little mass. 

Most important, this assumption is supported by direct obser- 
vations, which indicate that the central component is well-matched 
by a single stellar population (with the background galaxy - i.e. the 
violently relaxed component - being some older population), usu- 
ally better so than by mode ling it is as having resulted from star 
formation extended in time |Trager e t aL 2000; Nelan et al. 2005^; 
Thomas et alj|2005l; lOallazzi et alj|2005.. 2006 ; Lauer et al,.2005i ; 



Kuntschner et alj|2006l ; iMcDermid et alj(2006i) . The metaUicities 



and Q-enhancements of these regions require that the stars there 
formed rapidly, in short-lived events or in events closely spaced in 
time (later series of bursts, or independent well-spaced events, be- 
ing ruled out by the enrichment patterns). 

Regarding assumption (3); the observations argue that the 
system must have been (at least briefly) dominated by a strong 
gas inflow, that then formed stars more or less in situ. The short 
timescales of star formation indicated by stellar populations do not 
allow for a trickle of gas over an extended period of time. If as- 
sumption (1) applies, the corresponding star formation timescale 
required to match the abundances is close to that obtained if all 
the gas w e re in place at once, and formed stars according to the 
iKennicutlJ ( Il998h relation ~ i.e. the shortest allowable timescale. 
Moreover, in hydrodynamic simulations of galaxy-galaxy mergers 
or dissipational collapse, this is almost always the case - gas flows 
in on the local free-fall time. Ho wever, star format ion is inefficient 
on the the free-fall timescale (the lKennicutll ( Il998h relation imply- 
ing star formation efficiencies of ~ 3 — 10% per free-fall time), 
and so catches up over a short period following the initial gas-rich 
phase of inflow. In addition, the disky, rotational kinematics dis- 
tinctive of the burst relics requi re in situ formation from initially 
gas-rich configuratio ns (see e.g. lNaab et al ] |2006l; ICox et alj|2006l ; 
iHopkins et al.ll2009el) . 

Note that the assumptions outlined here clearly do not apply 
for stars not formed in bursts. Although observations suggest that 
the lKennicutll ( 1 19981) law applies (although there may be a surface 
density threshold for star formation), it is established, in contrast to 
the central, burst populations, that the extended, violently relaxed 
stars were formed through star formation histories more extended in 
time. Disks, of course, are fed through continuous accretion, rather 
than a single massive inflow; so both assumptions (2) and (3) break 
down. Moreover, the outer components of ellipticals are formed via 
violent relaxation - i.e. by definition from merging already-formed 
stars from other galaxies - and so have no reason to be a homo- 
geneous population. In fact, especially in massive ellipticals, it is 



cosmologically expected that such por tions of the galaxy are made 
up of debris from many small systems (Gallagher & Gstriker 1973 
iHernq uistetal. 19931; iDe Lucia etaU 2006; d e Luc i a & Blaizo3 
p07l ; iHopkins et aiTl2009j |il; iKhochfar & Silkl bood ; iNaab et all 
1 20091) . It is likely then, that these components - ~ 90% of the stars, 
in typical ~ L* systems - have star formation histories that are not 
trivially invertible. 

2.3 Methodology 

Together, these assumptions allow us to invert an observed 
Sbuist(^). If the relic burst mass is dominated by of order a sin- 
gle major event, and that event was indeed a dissipational, gas-rich 
and/or rapid event, then the gas density at the beginning of the burst 
(which we define for now, arbitrarily, as f = 0) should be given ap- 
proximately by 

Eg,,(/?,/=0) ^EburslW ■ (1) 

The iKennicut 3 l l 19981) law relates the star formation rate surface 
density - and hence the gas depletion rate - to the gas surface den- 
sity through 



'^-Eg,, ^ 1.5 X lO^^Mgyr^'kpc"^ 



r) (2) 



d/~"°' ' VMq pC-2. 

where the nor malization come s from IKennicut 3 l ll998h corrected 
for our adopted IChabrie3 ( |2003|) IMF, here chosen fixed for Milky- 
Way like systems (where the relation is best calibrated). The index 
iiK ~ 1.4 is suggested by observations of the local and intermediate- 
redshift Universe, and so we will adopt this as our fiducial value. 
However, a higher index of hk ~ 1 .7 has been suggested from ob- 
servations (albeit more u ncertain) of systems at high redshifts with 
extremely intense SFRs jBouche et alj|2007h ; we therefore com- 
pare our results to those obtained for such a higher iik- 

Note that since we are ultimately interested in the relic star 
formation rate and time-averaged gas depletion rate, stellar mass 
loss can be absorbed into the normalization of the SFR-surface 
density relation (in essence, adopting the instantaneous recycling 
approximation). But, in general, we fi nd that for reas onable recy- 
cling fractions from e.g. the models of iBruzual & Char iot (20oJ), 
the resulting differences are less than those between the choice of 
the slope iik. Likewise, the effects of long term stellar mass loss 
(although, post-burst, a typical IMF will yield just ~ 20% stellar 
mass loss) are generally smaller than uncertainties arising from our 
primary assumptions above. Fortunately, the uncertainties in both 
(e.g. the singularity of the burst versus degree of stellar mass loss) 
have opposing signs, so experimenting with reasonable variations 
in each yields little net difference in our predictions. 

With this relation, we can then integrate forward in time from 
the initial gas density in each annulus, computing new gas and stel- 
lar densities after some time interval, and then determining a new 
E, in that annulus. We show the results of this procedure for our il- 
lustrative example in Figure[T] At the initial time, the gas density is 
given by the surface density profile of the relic starburst component 
isolated from the observed profile (as seen in the Figure, it makes 
no difference if we take the empirical fit or starburst component 
derived frmo the best-fit simulation). As discussed above, we sub- 
tract out the non-starburst component. This leads to a correspond- 
ing SFR surface density. This depletes gas, leading to a lower gas 
surface density. The evolution is most rapid in the central regions, 
where the densities are highest, because the Kennicutt-Schmidt re- 
lation is super-linear in surface density. At several times, we plot 
the remaining gas and corresponding SFR surface densities, which 
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both decrease and become more extended as the dense, central re- 
gions become depleted. 

Integrating the SFR surface density over R at each time, we 
obtain the total SFR in the burst at f > 0, 



M. = 



(3) 



Moreover, if we assume that light, in the UV, radio, or IR, traces the 
star formation rate locally in some annulus, using e.g. the relation 



LiR ^ 1.14 X 10'"L(7 



■0 



(4) 



from 



KennicuttI ( Il998h (adjusted as appropriate for the IChabrieij 



that the burst would start exactly from pure gas ait — Q, as assumed 
here). Nevertheless, it is clear that these details lead to almost no 
deviation between the true SFR versus time and that inferred from 
our simple procedure, given the observed relic starburst/inner com- 
ponent of the system. This is in part because some of these effects 
are small. Also, several tend to offset one another, so that, on av- 
erage, the dominant physical considerations are the validity of an 
average Kennicutt-Schmidt law and the rapid, initially gas-rich na- 
ture of dissipational starbursts. 



3 RESULTS 



( l2003h IMF adopted here), we then we obtain the light distribution, 3.1 Typical Behavior: Some Examples 



and can evaluate e.g. the half-luminosity radius at each instant in 
time. ' 

Figure [T] shows the resulting total SFR as a function of time. 
Note that the procedure above, strictly speaking, defines only a one- 
sided SFR for f > - i.e. times later than some initial burst peak, 
M»(f > 0). Real starbursts, of course, have a rise as well as a fall 
after their peaks, and simulations (as well as simple dynamical con- 
siderations) indicate that - to lowest order - the rises and f alls are 
roughly symmetric jdi Matteo et al.ll2007l; ICox et allFlOOSh . If we 
assume such symmetry, it is trivial to convert our solution for the 
burst from f > to that for all times, as follows. 

Our modeling gives a one-sided solution for the SFR versus 
time at ? > 0; namely M*|one .sidcd(r > 0) = f{t > 0). Physically, 
our symmetry assumption means that half the time at each lumi- 
nosity is spent on each side of the peak. Quantitatively, for r < 
(pre-peak), the SFR M.(/ j f < 0) ^ f{-t/2), and for r > (post- 
peak), |f > 0) ^ f{+t/2). This is the symmetric or "two- 
sided" light curve. This choice guarantees that the total stellar mass 
formed in the burst is preserved, and enforces the symmetry con- 
straint. It is, therefore, appropriate for comparison to full time- 
dependent light-curves. However, we stress that for all the quan- 
tities derived in what follows, it makes no difference, whether we 
assume a one-sided or symmetric two-sided lightcurve. The dura- 
tion of the burst is the same. (If it is defined as some ?burst, we simply 
have to change our notation from referring to f = to ? = fburst, to 
t = — fburst/2 to f = +/burst/2; but the zero point of time is arbitrary 
in any case.) Moreover, for the purposes of luminosity functions 
and SFR densities, the relevant quantity is the time spent in each 
SFR interval - this is, by definition, identical whether or not we 
split that time across two sides of a peak. 

We compare the integrated SFR as a function of time from our 
illustrative example in Figure [T] inferred from the relic starburst 
profile, to the true simulation SFR as a function of time in the full 
hydrodynamic simulation that produces such a profile. The agree- 
ment is excellent, verifying our procedure. This is despite the fact 
that, in the full simulation, many of our assumptions are not true 
in detail. For example, the simulations obey a local (not global) 
Kennicutt-Schmidt law; thus clumping and instabilities can accel- 
erate star formation locally. The simulations also account for feed- 
back from stars, stellar winds, and black hole accretion, which can 
remove and recycle gas. The SFR is not precisely trivial or mono- 
tonic in time. And the gas inflows do not occur instantaneously (so 



' More sophisticated, luminosity and galaxy property-dependent conver- 
sions have been proposed, but they largely differ at low IR luminosities 
in e.g. extended disks. Since the systems of interest he re ai'e massive star- 
bursts , w e find that the alterna tive conversions from e.g. lBuat & Burgarellal 
i 19981) or ljonsson et al.l l200d) make no difference to our conclusions. 



Figure [2] illustrates the results of our starburst inversion pro cedure 
for a representative sub-sample of the lKormendv et al] l l2009l) Virgo 
ellipticals. We first show (top) the stellar mass surface density pro- 
files of the inner/dissipational/burst component in each object (not 
the total stellar mass densit y profile), as obtained using the two- 
component decomposition in lHopkins et al. l( l2009blB). over the dy- 
namic range covered by the observations from Kormendv et al.l 
( l2009h . In other words, we have already subtracted off the empiri- 
cally estimated non-burst stellar component. For each system, using 
the methodology described in § |2l we convert the inferred surface 
brightness profile to a time-dependent total star formation rate, as 
illustrated in Figure [T] We show the one-sided (t > 0) SFR versus 
time so obtained, in a logarithmic scale (this highlights the rele- 
vant power-law behavior). The total burst star formation rates decay 
from some initial maximum in a power-law like fashion. We find 
that good fits can be obtained to the time-dependent light curves 
with the functional form 



M.(f) = 



M,(0) 



M*(0) 



[l + {t/to)]f [l + (|f|/[fo/2])]'3 



(5) 



The first is a one-sided light curve appropriate if the burst occurs 
only for f > 0, and the second is a two-sided curve symmetric about 
t — (see § |2]l- The two functional forms shown are statistically 
equivalent for our purposes, and we therefore refer to them inter- 
changeably. Fitting such a curve to each profile in Figure |2] the 
best-fit curves are indistinguishable on the scale shown from the 
numerical calculations. We discuss the best-fit values below, but 
in general note that the systems tend to cluster around a power- 
law decay with /3 ~ 2.5, and clearly exhibit a characteristic decay 
timescale of ~ 0. 1 Gyr. 

Note that, with such a functional form, there are straightfor- 
ward relations among several parameters. The half-life of the star- 
burst ti i.e. the time for the starburst to decay to 1/2 of its maxi- 
mum luminosity/star formation rate, is simply 

tyl2 = {l'^" -\)t0. 



(6) 



Likewise, the peak SFR of the burst, M, (0), can be trivially re- 
lated (knowing the burst duration fi/2 and f3) to the total stellar 
mass formed in the burst, via the integral constraint that M^mA ~ 
/o°°M.(f)df, giving 



M.(0) = (/3-l) 



[(/3-l)(2 



1//3 



1)] 



tl/2 



(7) 



For convenience, we will use ti /2 as our timescale of interest below, 
as it has a straightforward interpretation independent of /3. More- 
over, for the /3 values of interest, the coefficient [(/? — 1) (2*'^'' — 1)] 
in Equation|7]is nearly constant at ~ 0.45 — 0.50. 

At each time, we can also estimate the effective radius of the 
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Figure 2. Typical results from recovering the burst star formation histo- 
ries of spheroids. Top: Surface stellar mass density profiles of the em- 
pirically fitted burst c omponent alone, for each observed galaxy in the 
iKormendv et al.U 2009h Virgo elliptical sample. Colors denote the total mass 
of each spheroid, from IO'Mq (violet) to ~ IO'^Mq (red). Middle: Im- 
plied burst SFR versus time since the peak of the starburst. The behavior is 
roughly self-similar and well-approximated by Equation |5] Bottom: Burst 
size scale (radius enclosing 1/2 the total star formation) versus time. A 
similar scaling applies (Equation[8). 



starburst region; we specifically define this as the half-luminosity 
radius, where we assume the surface luminosity density (in e.g. the 
infrared) is proportional to the gas surface density at each instant. 
We show this as well {bottom) in Figure [2] In analogy to the fitted 
total star formation rate, we can similarly describe this evolution by 
a power-law of the form 

RA{t)^RA{Q)[l + {t/to,r)f' . (8) 



The radius doubling time ti /2,r is trivially related to to.r by the same 
relation as Equation|6]above; i.e. ?i/2., ~ (2''''''" — l)fo.,. 



3.2 Properties of Bursts versus Mass and Scale 

3.2.1 Masses, Tiinescales, and Star Formation Rates 

Figures [3] illustrates these fits and a number of global proper- 
ties of the burs t s, usi ng o ur full sample of profiles taken from 
iKormendv et al. and lLauer et al Tl l2007h . We also present the 

relevant fits and scaling relations in Table [l] First, Figure [3] shows 
the total burst mass (i.e. total mass of the inner/dissipational com- 
ponent) as a function of total galaxy stellar mass {top left). These 
values, and independent tests of their validity, as well as compar- 
ison to models of spheroid format ion and obser ved disk gas frac- 
tions, Me_discussed extensively in iHopkins et al . (2008c, 2009bJ, 
l2008all20093) . To lowest order, a fraction ^ 10% is typically in the 
burst/dissipational component, i.e. Mbunsi '•^ 0. 1 M« . For som ewhat 
higher accuracy, we use the fit from lHopkins et aL I ( l2008ah , moti- 
vated by comparison with disk gas fractions that constitute burst 
progenitors. 



(9) 



[1-f (m./io'-'^Mq)"*] ■ 

This is also shown in Figure[3] along with the simpler linear fit. 

Next, we show the burst half-life ti^2 {bottom left), as a func- 
tion of burst mass (the comparison versus total stellar mass is sim- 
ilar, but other possible correlations are less significant). To lowest 
order this is mass-independent, with a value ~ 10** yr - depending 
on the definition, this equates to a total observable burst lifetime of 
a couple 10* yr with an inherent factor 2 — 3 scatter. Fitting to a 
power law, we find some (weak) mass dependence. 



'1/2 ' 



1.1 X 10 yr 



(10) 



: 0.2 dex about this me- 



with a mass-independent scatter of at^ 
dian relation. 

This agrees well with expectations from e.g. n umerical sim- 
ulations and extrapolation of the Kenniciit^ i ll 9981) law, that the 
burst duration should reflect the dynamical times in the central re- 
gions of galaxies; as e.g. a disk galaxy dynamical time oc Re/Vc 
is only weakly mass-dependent (oc yvft"-"^""-"^)^ jf yy e compare 



Re OC A^^' and Vc (x M°-^^ from 'Courteau et alj|2007t) (see als o 



iBell & de Jon3 1200 ll ; iMcGauglj i2005i ; .Avila-Reese et al.l l2008h . 
Moreover, the scatter in the dynamical times expected from these 
correlations is similar to that in Figure[3] 

We next show the slope of the power-law decay, /3 {bottom 
right). There is no obvious correlation with burst mass or any 
other parameter; however, the values cluster in a reasonably nar- 
row range, /3 2 — 3, with a best-fit median of /3 ~ 2.38. Again, 
this matches expectations; given some Schmidt-law star formation 
relation, the slope /3 follows from the profile shape of the relic 
burst component. These pro file shapes, and their physical origins, 
are discussed extensively in iHopkins et al. Il l2008cl , l2009bl) ; there, 
the authors show that they can be reasonably well-parameterized 
as Sersic functions with relatively low Sersic indices n, --^ 1 — 3, 
and that this follows naturally from the fact that they form via gas 
inflow and dissipation. Given such a Sersic-index profile, although 
there are no trivial analytic solutions for (3, a simple numerical cal- 
culation shows that /3 ~ 2 — 3 is the natural, corresponding expec- 
tation. 

We next consider the peak burst star formation rate, M*(0) 
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Figure 3. Observationally inferred properties of splieroid starbursts as a function of mass, from the relic starbursts. Top Left: Mass of tlie starburst component 
of spheroids as a function of total spheroid stellar mass. Detailed d eterminations and tests of these component de compositions are pr esented in H09b,e. We 
show results for the Virgo galaxy sample o f lKormendv etal] i2009l red stars) and local massive spheroid sample of lLauer et alj lioovl violet circles). Dashed 
blue lines show the range of best-fit scalings to the observed systems (given in Table[T). Top Right: Peak SFR of the starburst, as a function of total starburst 
mass. Bottom Left: Duration of the starburst (time from peak to half-peak SFR). Bottom Right: Best-fit power-law slope to the decay of the SFR versus time 
(Equation[5j. 



Table 1. Starburst Scaling Relations 



.V 




X 


a 


b 




log(Mb„,,,/IO"'MQ) 


\o% 


(M,/10"Mo) 


-0.05±0.12 


0.78±0.I6 


0.41 


log(M.[0]/MQyr-l) 


\o% 


(Mb^^t/IO'OMQ) 


1.70 ±0.05 


0.92 ± 0.04 


0.22 


log ((i/2/Gyr) 


\o% 


(Mb,,,t/10'0MQ) 


-0.96 ±0.04 


0.08 ± 0.04 


0.20 


log(M.[0]) 


log 




-0.24 ±0.01 


1.00 ±0.01 


0.04 




\o% 


{Mb„„,/10"'MQ) 


2.38±0.12 


0.09 ±0.11 


0.35 


log (Rrdic/kpc) 


\o% 


{Mbu^t/lOlOMg) 


-0.43 ±0.07 


0.58 ±0.07 


0.29 


log(Ri/2/kpc) 


\o% 


(Mbu^t/lQlOMQ) 


-0.57 ±0.07 


0.60 ±0.06 


0.26 


log (fi/2.r/Gyr) 


\o% 


(Mb^^t/lOlOMg) 


-1.07 ±0.08 


0.10 ±0.08 


0.28 




log 


(Mb.^.t/lOlOMQ) 


0.45 ±0.01 


-0.035 ±0.017 


0.11 


log (Burst Age / Gyr) 


\o% 


(M./IQI'Mq) 


0.95 ±0.11 


0.15 ±0.04 


0.12 


Burst Age/Gyr 


log (M*/1O"M0) 


9.5 ±0.4 


3.0 ±0.6 


2.1 



For the parameters x and v, the best-fit linear scaling of y(.v) is presented, with the assumed 
functional form y = a + bx. EiTors on the fitted parameters a (normalization) and b (slope) 
are given, along with the intrinsic scatter in y about this correlation, a. 
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{top right). Immediately, it is clear that this is very tightly corre- 
lated with the total burst mass Mbmst- We compare the correlation 
expected from Equation [7] assuming a typical /3 ~ 2.4 such that 
(/?- 1) (2'/" -l)^ 0.5, i.e. M, (0) f» O.SMbursi Ai/2. We show this 
both for an assumed constant /[/t = 5 x lO'yr, and for /i/t given 
by the weakly mass-dependent fit in Equation[TO] Both cases agree 
very well with the observed trend, with a small, mass-independent 
scatter of a^^^ ~ 0. 15 — 0.2 dex. 

Worth noting is that, within massive observed systems, peak 
burst SFRs of > 1000 MQyr^' are obtained, but there is a large 
dynamic range ~ low mass systems can have "bursts" with implied 
peak SFRs as low as '--^ 0. 1 — 1 Mq yr^ ' . Also, since there is rel- 
atively little variation in ti/2 at fixed mass, the scatter is, as men- 
tioned above, small. Producing a very high-SFR burst therefore re- 
quires an extremely high-gas-mass system; the case where a lower- 
mass system might compress its gas much more efficiently, leading 
to a shorter-lived but arbitrarily highly-peaked burst appears to be 
rare. 

Finally, we note that their is some (factor ~ 2) observational 
uncertainty in the exact normalization of the correlation between 
SFR and gas surface densities. Adopting a different IMF will 
likewi se systematically shift the inferred relation from iKennicuttI 
( Il998l) . However, so long as the change in Equation [2] is purely 
in normalization, it has no effect on the shape of the inferred 
lightcurves, and translates directly to corresponding normalization 
changes in the resulting SFR and timescales in Figures l2P1 and Ta- 
ble [T] Specifically, if the normalization in Equation |2] is multiplied 
by a factor -q, then the implied peak SFR M, [0] increases by the 
same factor i], and the burst timescale fi/2 (and ti/2.r) decrease by 
a factor rj. All other parameters in Table[T]are unchanged. 

3.2.2 Spatial Sizes of Bursts 

Figure |4] continues our analysis, plotting parameters as in Figure [3] 
but for the burst spatial sizes. First, we consider the effective (pro- 
jected half-stellar mass) radius of the burst remnant (top left), i.e. 
the radius of all stars after the burst has completed, namely of the 
relic profile as shown in Figure |2] (fop left). We consider this as a 
function of burst mass, and find a correlation of the form 

...^0.37kpc(^)"-'^ (11) 

with scatter o-r^^^.^ ~ 0.3 dex. As a p roperty of the total mass distri - 
bution, it is s tudied in more detail i n lHopkins et ah Ii l2008d . l2009bh . 
There, and in iHopkins et a"i]j2009il) , it is shown that such sizes arise 
naturally from dissipational inflows, as a consequence of where the 
inflows stall and become self-gravitating. This is the reason that the 
power-law scaling in Equation ll Ilreflects the obs erved scaling of 
spheroid sizes with mass (see e.g. lShen et al .l2003h . and is a conse- 
quence of simple dynamics, independent of the star formation law. 

Next, we consider ^bmsi {top right), the effective {half-light, 
assuming light traces local SFR) radius of the burst itself during 
the starburst episode (as shown in Figure |2j bottom). For conve- 
nience, and as a reference point, we define the plotted quantity as 
/?bursi {ti 12 ) , i.e. the size at time / = /1/2 . This is a convenient param- 
eterization and also has the advantage of being more representative 
of what might be observed, since the duty cycle at peak is vanish- 
ingly small. In any case, there is a similar correlation, with 

/ M, , \ 0.60 

R,I2 = Ru.4h/2) ^ 0.27kpc (^qT^^^) (12) 
with o"R„,j^. ~ 0.3 dex scatter. 



As expected, the size Rbmn near the peak of star formation is 
much smaller (by a typical factor ~a few) than the size scale of 
the relic stars. This owes to two factors. First, near-peak, it is the 
highest-density material just at the galaxy center that dominates the 
SFR, while almost none of the more extended material is contribut- 
ing at a significant level - hence the size of the burst itself will still 
grow by a large factor with time (see Figure |2j. Second, this size 
is weighted by S FR, i.e. by gas de nsity to some super-linear power 
according to the IKennicuttI ( Il998h relation. If Egas follows a Sersic 
law with some index and effective radius Re, and Esf oc E"fs, 
then it is trivial to show that the half-SFR radius will simply be 
Re X iij^"' . For lis — 4 and iik = 1.5, this gives a factor of almost 
5 smaller Re for the SFR distribution (for = 2, more typical of 
the inner components of the systems here, this is a more moderate 
factor ~ 2). 

Nevertheless, the spatial sizes can still reach several kpc. In 
particular, for the same systems that reach > 1000 MQyr^' peak 
SFRs, the effective radii are also large, ~ 1 — 10 kpc. Moreover, we 
refer here only to the strict half-SFR radius; estimating an observed 
size in detail requires radiative transfer modeling. Analysis of simu- 
lations in IWuvts et al. ( 2009b , a ), and similar analys is of some well- 
studied local starburst galaxies l lLaine et alj2003bl) . suggests that - 
especially at the peak of activity - the IR or radio size can be biased 
to factors of ~ 2 larger than the true half-SFR radius owing to ob- 
scuration in the high-column-density nucleus and reprocessing of 
the emitted light. 

Next, we consider the burst size-doubling timescale, fi/2.;. 
Unsurprisingly, the values are very similar to the SFR decay times 
ti 12, with similar weak dependence on mass. We find a best-fit me- 
dian/1/2,,. 

V2.-8.5xl0yr(^^) , (13) 

with a mass-independent scatter of (t,|^, . ~ 0.3 dex. 

We also consider the slope /3, , of the starburst size versus time. 
Here, there may be a weak correlation with mass at the high-mass 
end, but formally this is only marginally significant (~ 2a) and 
only if we restrict ourselves to Mbmst ~ lO'" — 10'"Mq. We there- 
fore hesitate to quote a correlation; but do note the relatively ro- 
bust median value of /3r ~ 0.4 — 0.5, with scatter tj/j,. ~ 0. 1. Again, 
the se naturally follo w from the combination of profile shapes and 
the lKennicut T h998l) relation. 

3.2.3 Effects of a Steeper Kennicutt-Schmidt Relation 

Thus far, we have adopted the Kennicutt-Schmidt relation as cal- 
ibrated at low redshifts (Equation with an index iik = 1.4 
(E, oc EjJs)- However, some recent observations of high-redshift, 
massively star-forming systems have suggested that the SFRs of 
such systems exceed the predictions with th is index and favor a 
somewhat steeper slope of iik = 1.7 ±0.1 jBouche et al.l 120071) . 
Specifically, they argue for a best-fit relation 

E* ^9.3 x 10"^M(7)yr"'kpc"^ ( ^ V'^ . (14) 

VMqpc"-/ 

We caution that the high-redshift observations - in particular the in- 
ferences of both gas masses and total bolometric luminosities and 
hence SFRs - remain uncertain. Nevertheless, this is an important 
source of error, and generally larger than many of the other sys- 
tematic uncertainties (those being at the factor of ~ 2 level) in our 
modeling. Moreover, it can be particularly important in the most 
extreme systems, such as e.g. sub-millimeter galaxies, with high 



© 0000 RAS, MNRAS 000, 000-000 



10 Hopkins et al. 



Best-Fit Trend 




0.01 



7 8 9 10 11 
log( Burst Mass M^..^, ) [ M®] 



0.01 



7 8 9 10 11 
log( Burst Mass M^.^s, ) [ M®] 



10 




0Q 



O 



o Lauer et al. Sample 
i!VKormendy et al. Sample 



8 9 10 11 
log( Burst Mass M^,,,,^ ) [ M®] 



12 




1.0 

CD 

8- 0.8 

CO 



0.4 



o 
O 

i) 

N 

CO 

^ 0.2 

CO 



0.0 



o 



o o 



o 




12 



8 9 10 11 
log( Burst Mass M^..^, ) [ M®] 



12 



Figure 4. As Figure [3]- observationally inferred properties of spheroid starbursts as a function of mass, from the relic starburst stars - but here giving the 
spatial size distribution of starbursts. Top Left: Size (projected half-stellar-mass radius) of the burst relics (generally a factor ~ 0.1 of the total galaxy size Rg), 
directly from the observed galaxy properties. Top Right: Size of the starburst at time t = 1^/2 (when the SFR is half-peak). Size here is defined as half-SFR or 
half-hght (assuming light traces SFR) radius; hence is more compact than the relic size. Bottom Left: Size-doubling time of the starburst (time for e.g. relative 
gas exhaustion in central, high density regions to leave only star formation at larger radii, giving a larger apparent radius). This is simply related to the SFR 
decay timescale in Figure|3] Bottom Right: Best-fit power-law slope to the increase of starburst scale radius versus time (Equation[8j. 



gas densities. We therefore re-consider our previous comparisons, 
but instead adopt this modified (steeper) index. 

Table|2]shows the resulting revised fits, for quantities that are 
affected by this (e.g. the SFR and size versus time). The resuhs 
are, in all cases, qualitatively very similar to those in Figures [2P1 
We therefore show only the results of fitting the correlations, and 
present a simple illustration of how the inferred burst properties 
change with iik in Figure |5] Specifically, we repeat our analysis 
of the observed Virgo systems from Figure [2] but include both the 
previous results (with iik — 1 .4) and the steep-index results (with 
hk = 1.7). For clarity, we show only a random subset of the ob- 
served systems. There are significant, systematic offsets between 
the implied star formation histories. 

As expected, with a steeper index - i.e. larger SFR at high den- 
sities, the peak SFR of bursts is systematically higher. The differ- 
ence is large - a factor of ~ 5 higher peak SFR at fixed burst mass 
- giving (M.(0)) ^ AIOMqyi-' {Mb,r.t/W'°MQf^'^- in other 
words, this implies that most systems with Mbmst > 3 x IO'^Mq 
will exceed 1000 Mq yr^ ' SFRs at the peak of their starburst phase. 
However, with a higher SFR comes faster gas exhaustion, so the 
half-life of this peak, fi/2. is systematically shorter, by a factor of 
~ 6 - i.e. typical exhaustion times ti/2 ~ lO'yr. The decay slope 



/3 is also systematically different, now /3 ~ 1.7 (as opposed to 
2.4 before) - this shallower power-law decay may appear counter- 
intuitive, but, from the form of Equation|5] relates to the fact that the 
SFR at r ^ ti/2 is more peaked with iik ~ 1.7, and then at f 3> fi/2 
(when the gas density is lower) more extended (because the steeper 
index yields lower SFRs at low surface densities). 

In practice, if one defines a starburst "duration" by time spent 
above some fixed, high SFR, this change in the Kennicutt-Schmidt 
slope does not actually have as dramatic an effect as implied by 
the change in fi/2. At f > 0.1 Gyr, in fact, there is relatively little 
difference in the SFR versus time, starting from the same observed 
relic starburst density profile (i.e. for the same relic object). Most 
of the difference comes at f < 0. 1 Gyr. This is clear in Figure[5] 

3.3 Inferring the Burst Star Formation History 

3.3.1 Methodology: The Ages of Bursts Today 

Having quantified how burst properties scale with mass, and how 
these are connected to overall properties of their host galaxies, we 
require only a simple extension of our analysis to infer the cosmic 
history of starbursts. 

Figure [6] shows the observed ages, at z = 0, of the domi- 
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Table 2. Starburst Scaling Relations with Alternative Kennicutt-Schmidt Relation 
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0.11 ±0.09 


0.30 


log 


(M*[0]) 
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log 
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-2.53 ±0.23 


0.06 ±0.21 


0.41 




log(A/bur.„/lO"'M0) 


0.33 ±0.01 


-0.02 ±0.01 


0.08 



As Table[T] but adopting a Kennicutt-Schmidt star formation versus gas surface density 
relation with a steeper index (S* oc Sg^j), as suggested by obsei-vations of high-redshift, 
high-SFR systems. 
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Figure 5. As Figure |2] {middle), but comparing a subset of the inferred 
star formation histories from the observed systems with different assumed 
Kennicutt-Schmidt law indices. Solid lines, as Figure |2] assume the z = 
index ha: ~ 1 .4. Dashed lines assume a steeper relation suggested by some 
high-redshift observations, riK ~ 17. The latter leads to more sharply 
peaked bursts with larger peak SFRs in early (t < 10** yr) stages. The star 
formation histories after ~ 0.1 Gyr are similar, once gas surface densities 
have decayed to more moderate levels. Table|2]shows the results from re- 
fitting the correlations in Figures |3j4] with this higher nx = 17. 



nant starburst in z = ellipticals/spheroids as a function of their 
total ste llar mass. We compile observations from a number of 
sources dXrager et al.l20o"ol:ICaldwell et alJ2003l:lNelan et al.l200p : 



iThomas et al. 2005k iGallazzi et all l200a iBemardi et al.l I2OO6I1 
these involve various observational methodologies and selection 
functions, and together should plausibly represent the systematic 
uncertainties in this process. Note that for some of these sam- 
ples, ages are quantified by environment; since the mass den- 
sity of spheroids at the masses plotted is dominated by those in 
field environments, we show these values. Other observational es- 
timates are within the scatter defined by these observations (e 



Ijeirgensen et al.ll999l : lKuntschner et aLbOOlhlGaTlazzi et alj2005l 
For each sample, the median age is well-determined and agrees rea- 
sonably well with other estimates. A good approximation to the 
median ages is given by 

M* \0-i3 



fiookback.biiisi ~ 8.9Gyr 



10" M, 







(15) 



We also plot the 1 a scatter in ages at each mass determined from 
each sample; these too agree well. We find that this is reasonably 
parameterized as a mass-independent, Gaussian ±2 Gyr scatter. 




□ Trager et al. 2000 
A Caldwell et al. 2003 
oNelan etal. 2005 
V Thomas et al. 2005 , 
■^Bernard! et al. 2006 Z 
o Gallazzi et al. 2006 
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Figure 6. Observed average single stellar population age (at z = 0) of 
spheroid nuclear/most recent burst populations, as a function of total 
spheroid stellar mass. Points show the median age at each mass from various 
determinations; error bars show the scatter (not the much smaller en'or in 
the mean). Observations are compiled from various sources (labeled). Solid 
black lines show the best-fit log-log and log-linear trends for the median 
values (Equation[T5]and Table[T). 



Various analyses, for example observationally analyzing mock 
SEDs constructed from hydrodynamic galaxy formation simu- 
lations or semi-analy tic models of galaxy form ation (see e.g. 
IWuvts et al.i ,2009b, a.; IXrager & Somervilld I2OO9I , and references 
therein) have shown that, for typical spheroids that combine stars 
formed over an extended (potentially still short in absolute terms) 
star formation history in pre-merger disks (assembled dissipation- 
lessly in mergers) with stars formed in bursts via gas dissipation in 
mergers, these observed ages do reflect the cosmic time at which 
those bursts occurred. Since a number of constraints and theoret- 
ical models show that the time of last major gas-rich merger and 
time when star formation in "quenched" systems must shut down 
are tightly coupled (even though there may not necessarily be a 
causal relationship between the two), this can also be thought of 
as dating the time of the last burst of star formation that con- 

^2008 

Haiman et a l. 2007; Hopkins et al. 2007, 2006a; 



20061 ; iHopkinset al.ll2008bl ; IShankar et al.ll2009l) . Moreover, most 



of these ages are derived specifically from the central stellar popu- 
lations, those that are clearly formed dissipationally; thus they rep- 
resent a clear constraint on when this burst occurred independent of 
more complex assembly and star formation histories for the more 
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extended components. Finally, where available, detailed, resolved 
measurements of stellar population gradients show that these cen- 
tral populations do represent a distinct component, in agreement 
with dissipational models, and that their properties favor a single 
burst, rather than an extended star formation hist ory or dissipation- 
less assembly of many distinct sub-components jSchweizedl 19961: 
Titus et al.' 1997:'Sc hweizer & Seitzeilll998ll2007l : lReichardt eraU 
2Q01. ; .Mich ard 20061). 



Given these comparisons, it seems reasonable to take the age 
distributions in Figure[6]as a first approximation to the times when 
these starbursts occurred. Together with the analytic fits to the dis- 
tributions of burst properties versus mass, this enables us to recon- 
struct the burst history of the Universe. Specifically, we begin with 
the mass functi on of bulges/sphe roids at z = 0, here adopting the 
determination in iBell etal.l 120031) . We ignore disksAate-type galax- 
ies and very low-mass bulges, as these contribute negligibly to the 
quantities of interest (they have very low burst masses, hence low 
peak SFRs, where we will show the population is dominated by 
non-burst star formation).^ At a given mass, we know the num- 
ber density of systems today, and the distribution of burst masses 
Mbuist(M,,tot) (from e.g. Equation[9l with appropriate scatter); i.e. 
we know the (co-moving) total number density of bursts with var- 
ious masses that must have occurred before z = to account for 
those systems. We also know the distribution of cosmic look-back 
times when those bursts happened, from Equation[T5]above (again, 
with appropriate scatter). This allows us to construct a Monte-Carlo 
population of bursts, each with some total mass and taking place 
(defined here as peaking) at some cosmic lookback time. Note that 
this is independent of the state of the remainder of the galaxy at 
this time - so is not affected if e. g. the outer portions of ellipticals 
are assembled much later (s ee e.g. lNaab et al.l20()3 : lBezanson et al.l 
l2009l : lHoDkins et al.ll2009dllil) . 

We can then assign a full lightcurve to each burst, using the 
relations in § 13.21 and Table [T] namely, knowing the burst mass, 
a duration ti/2 (Equation 11 Ot and power law decay slope (median 
« 2.4; again with scatter). Together these define the correspond- 
ing peak SFR M, (0). The time f = in the lightcurve is simply 
the time of the burst, as determined above, and it evolves forward 
from this time (the exact choice makes little difference, since the 
burst durations are small compared to the Hubble time). For conve- 
nience, we truncate each burst at a time ~ 2O/1/2, but by this point 
the luminosity has decayed sufficiently that the burst is negligible. 
We then make mock observation of the Monte Carlo population at 
any redshift z, and construct the luminosity function or luminosity 
density of bursts. Of course, we are really predicting a distribu- 
tion of star formation rates; for comparison with observations, we 
convert these to total (8 — lOOO^m) infrared luminosities with the 
empirically-calibrated conversion factor in Equation|4] ' 



^ We have experimented with altamative determinations of the bulge- 
dominated galaxy m ass f unction at z = 0, including those presented in 



iBernardi et all i2009l) and Ivika et all ilOOfk . We have also attempted to 
constr uct the mass function of bulges specifically, following iDriver et alj 
i2007t) or by adopting the morphology-separated mass f unctions from 
iKochanek et alj )200ll) with a type-dependent B/T from iBalcells et alj 
)2007bl) . We find these make little difference to any of our conclusions, as 
the resulting differences lie primarily in identification of low mass bulges or 
bulges in late-type galaxies which are negligible in the total IR luminosity 
density. 

^ Note that, in comparing to observed starburst light curves or luminosity 
functions, the average M» / L will depend systematically on the IMF. How- 
ever, the inferred stellar surface mass density on which we base our analysis 



We can also trivially repeat this procedure for a different bolo- 
metric correction from M, to Lir (which will systematically shift 
the predicted luminosities by some factor ~ 2). Likewise, we can 
repeat our calculation using the fitted scaling relations appropri- 
ate for a steeper Kennicutt-Schmidt law with index «a- ~ 1.7 rather 
than iiK ~ 1.4, given in Table|2] 



3.3.2 The Luminosity Function of Bursts 

Figure|7]shows the resulting burst luminosity functions at a variety 
of redshifts from z = — 4. We compare with observations com- 
piled from a number of sources, available from z '--^ — 3. Note that 
all of these are corrected to a total IR luminosity from observations 
in som e band; we adopt the corrections compiled iniValiante et al] 
( l2009l) . but emphasize that some caution, and at least a systematic 
factor ~ 2 uncertainty in Lir, should be considered in estimates 
from most if not all observed wavelengths. 

At low luminosities, the predicted LFs are well below 
those observed, but they grow rapidly in importance, and are 
roughly consistent with the observations, at the high-luminosity 
end. This is expected - it is well-established that star forma- 
tion at relatively low rates is dominated by quiescent star for- 
mation in normal (e.g. non-merging) galaxies; i.e. distributed 
star formation in disks , rathe r than dissipational s tarbursts (see 
e.g. ISanders & Mirabell Il996l; iTacconi et al.1 12OO8I; iNoeske et al.l 



l2007b"a': 'Bell et al."2005'; 'Jogee et alj |2009l ; iRobaina et al.l 120091 
Veilleux et al. 2009; Sajina et al. 2007). At low redshifts, for exam- 
ple, nuclear starbursts (typically merger-induced) are negligible at 
luminosities < 10' ' ; but at the highest luminos ities, > 10'"Lq, 
they become dominant jSanders & Mirabell 1996h . At higher red- 
shifts, the luminosity function of bursts increases rapidly, as does 
the global luminosity function. This again is expected, as increas- 
ing gas fractions and specific star formation rates lead all sys- 
tems (mergers and quiescent disks) to higher SFRs at fixed mass. 
Thus, at all redshifts the (relatively) low-luminosity population re- 
mains non-burst dominated, and the threshold (in terms of Lir) for 
burst domination moves up to higher luminosities (Lir > lO'^ Lq at 
z ~ 1, and Lir > lO'^ Lq at z ~ 2). In terms of space densities, the 
burst population becomes important at ~ lO^^Mpc^"* log"' Lir, 
with much weaker redshift dependence. This is similar to the be- 
havior pred icted in models of merger-induced star formation bursts 
l lHopkins et al. 20 09k). 

We show predictions for both the cases of the low-redshift 
Kennicutt-Schmidt law slope of nK ~ 1.4, and for the suggested 
steeper value of «a- ~ 1.7 from high-redshift observations. Un- 
surprisingly, assuming a steeper relationship implies more intense 
SFRs at the peak of activity, so the predicted burst LF extends 
to higher luminosities. However, the shorter gas exhaustion time 
means that there is a lower number density at low/intermediate lu- 
minosities. The difference is clearly non-trivial, given the rarity of 
the highest-luminosity sources. The results in the steep-index case 
are comparable, however, to those obtained from the shallow-index 
case if we were to assume the scatter in the relevant quantities is 
systematically larger by ~ 0. 1 dex - in other words, this continues 
the same trends, and the effects may be comparable in magnitude 
to other uncertainties discussed here. 

We also provide simple fits to the inferred luminosity function 



is determined from observed luminosity profiles, with the same systematic 
dependence on M, /L. Thus, so long as the IMF does not evolve with red- 
shift, variations between typical choices make little difference. 
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Figure 7. Total (8 — 1000 /im) IR luminosity functions as a function of redshift. We show the contribution from bursts, using the methodology presented here 
(black lines), assuming a Kennicutt-Schmidt index of tin = 1 .4 (solid) or 1 .7 (dashed). Dotted red hne shows the result neglecting the scatter in the correlations 
shown in Figures |3P1 including the obsei'ved dispe rsion is import ant to expl ain rare, extreme stai'bursts. Points show observed estimates fromJSaundersetaT 



199C , magenta stars).ISoifer & Neugebauei] )l99lL blue triangles). lYun et aU ("200 1'. black filled circles), Huang et al.'l'2007', magenta squares), |Le Floc'h et aL 
2005 , violet diamonds), Caputi et al. 1 20071 orange inver ted triangles). jMagnelli et al. (200^ black *'s) Babbedge et al. (2006, red pentagons). Chapman et al 
20051 dark green +'s), and lPerez-Gonzli^^ m blue open circles). Shaded green rang e shows the prediction (with systematic uncertainty) for the IR 



LF of merger-induced bursts from the semi-empirical models and hydrodynamic simulations in lHopkins et al . ( 2009k JS. The burst contribution dominates the 
bright end of the IR LF, agreeing well with predicted merger-induced bursts and other constraints, but is a small fraction of the typical ~ L, activity at any 
redshift. The luminosity threshold above which bursts are important increases with redshift, along with the entire LF. Fits to the burst LF inferred here are 
provided in Table[3] 



of bursts, as a function of redshift, for both the shallow and steep 
Kennicutt-Schmidt slopes. Following standard convention, we fit 
the luminosity function with a double power law form, 

= 'tl (16) 

dlogL {L/UY + {L/U)^ 

where the parameters 0* (normalization), L, (break luminosity), a 
(faint-end slope, i.e. $ oc L^™ for L<^L,), and 7 (bright-end slope, 
i.e. $ oc L^^ for L 2> L^,) depend on redshift, with that dependence 
conveniently approximated as 



logL, 




log(^* 




a 


— ao + a ^ + a 4 


7 


= 70 + 7 ?+7 ? 




= log(l+z) 



(Note that log here and throughout refers to log|Q.) We perform this 
fit using our results only up to redshift z = 4, as the uncertainties 
grow rapidly at higher redshift. The fits should be considered with 
caution at higher redshifts. We provide the best-fit parameters in 
Table[3] along with their uncertainties. 

We note briefly that, especially because bursts are important 
at the high-luminosity end of the total LF, incorporating the scat- 
ter in the relevant relationships in § l3.2l is critical to the predicted 



abundance of IR-bright systems. In Figure |7] we compare the pre- 
diction if we were to ignore all scatter - i.e. simply construct our 
Monte Carlo population using just the median values of all param- 
eters and their correlations. As expected, this cuts off much more 
quickly at the high-luminosity end, reflecting the rapid exponential 
cutoff in the o bserved number of hi gh-mass galaxies. 

Recently, iHopkins et al.l ( l2009kl) presented predictions from 
galaxy formation models for the contribution to the SFR and IR lu- 
minosity distributions from normal, quiescent star-forming galax- 
ies, from merger-induced bursts of star formation, and from AGN. 
The models used a halo-occupation based approach to populate 
galaxies at each redshift (i.e. simply beginning with observed 
galaxy stellar mass functions and gas masses, with merger rates 
determined from evolving this forward in agreement with observed 
merger fractions), then mapping each population to suites of high- 
resolution simulations, to predict the distribution of SFRs in vari- 
ous systems. We compare their predictions for the merger-induced 
starburst population to our inferred burst SFR and IR luminosity 
distributions in Figure [7]'* Given the systematic uncertainties they 
quote (shown as the shaded range in the figure), and ours here, the 

^ The merge r rates determined in this model ai'e presented in 
iHopkins et iO )20091|) . A "merger rate calculator" script to give the 
merger rate as a function of galaxy mass, mass ratio, and gas frac- 
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agreement is reasonable. The two diverge at z 2> 3, but this is where 
the uncertainties in both the modeling and our empirical inferences 
become large. The faint-end extrapolations are also different, but 
in neither case are these well-constrained (in a systematic sense, at 
very low SFR/late times, the designation as "burst" is somewhat ar- 
bitrary). In general, the agreement seen supports our interpretation 
of the burst components of galaxies, and favors a possible merger 
origin for these bursts. More important, the burst history inferred 
here should not correspond to the normal or quiescent star forming 
population. 



3.3.3 The Problem of Bright SMGs 

We also note that, at z ~ 2 — 3, our reconstruction appears to some- 
what under-predict the abundance of the most luminous starburst 
systems relative to observations. However, at this redshift, the only 
constraints at the high luminosities of interest com e from the sub- 
millimeter populations observed in IChapman et"ZI feOOSi) . This is 
a relatively small sample, with a number of difficult completeness 
corrections involved in estimating the number density, and non- 
trivial cosmic variance given the sample selection. Recently, analy- 
ses of the number counts of similarly bright sources in much larger 
IR surveys have suggested that the average counts may actually be 
much lower - a factor ~ 5 difference ( with higher es timates owing 
to substantial cosmic variance; see e.g. lAustermann et al. 2009). 

Perhaps more important, given that the number densities in 
this regime fall rapidly with Lir, a small error in the bolomet- 
ric correction translates to order-of-magnitude differences in num- 
ber density at fixed L. And indeed, such a conversion from the 
sub-millimeter to total IR flux is highly uncertain, depending on 
quantities that are not well-known for these populations, such as 
the dust temperature, and subject to large object-to-object varia- 
tion in the local Universe. It should also be stressed that the ob- 
served high-redshift points in Figure |7] have not been connected for 
possible AGN contributions to the IR luminosity, although (from 
indirect constraints) it has been argued that AGN are unlikely 
to contribute more than ~ 30% of the bolometric luminosity in 
these systems i Menendez-Delmestre et alj2009l ; ICasev et alj2009l ; 
iBussmann et alj2009l) . For these reasons, we consider the compar- 
ison at these redshifts and luminosities to be largely qualitative. 
A more detailed comparison would require full modeling from e.g. 
high-resolution simulations that include spatially dependent, multi- 
phase star formation, metallicity, gas, and dust distributions that 
can forward-model the full SE P of such bur s ts (for some prelim- 
inary such results, we re fer tolLi et al.ll2008(; Ijonsson et al.l 1 20091 ; 



iNaravanan et alj2009bl lal; lYounger et alj2009l) . 

But the comparison in Figure |7] emphasizes an important 
point; ;/ the number density and bolometric correctio ns of such 
high-l uminosity systems are correct, then either the iKennicutj 
( Il998r) law must break down severely at high redshifts (giving 
much higher star formatio n rates than implied by e.g. observations 
of high-density systems in lBouche et alj2007l . at fixed surface den- 
sity), or something must be fundamentally wrong in our conver- 
sion between mass and light, either at low redshifts (i.e. some dra- 
matic errors in measurements of the surface brightness and stellar 
masses of local ellipticals, which appears unlikely) or at high red- 
shifts (i.e. the results of a strongly time-dependent stellar IMF, or 



tion, which determines these luminosity functions, is publicly available at 
|http : / /www .cfa.harvard. edu/~phopkins/ Site/mergercalc 



large heavily-obscured AGN contributions to the high-z luminosi- 
ties). 

Taken at face value, the number density of high-z bursts in- 
ferred, coupled with their implied star formation rates, would imply 
far too much high-density material in massive systems today. All 
of our analyses of burst star formation rates effectively set an upper 
limit on the burst number density at high-L. Consider, for example, 
the consequences of relaxing our assumptions. If starbursts were 
split into several separate events, with lower initial burst rates, then 
the duty cycle would go up, but the absolute L would go down (and 
in this portion of the luminosity function, the net result would be a 
severe decrease in the high-LiR prediction). If systems were not ini- 
tially gas-dominated or 100% gas, but maintained some lower gas 
fraction in quasi steady-state for the time needed to build up their 
densities, then again the peak luminosities decrease severely. If the 
present-day densities are the result of assembling many systems, 
this is the same as breaking up bursts, and leads to a lower pre- 
diction. The bursts cannot be more co ncentrated i n tim e (yielding 
higher peak SFR) without violating the lKennicuttI ( Il998l) law; even 
doing so, they become more luminous but shorter, and it requires 
an order-of-magnitude change to the star formation relation before 
this yields a match. 

Ultimately, a number of other obvious comparisons make this 
clear. Yielding the implied SFRs of the observe d bright SMG sys- 
tems (taken at face value) without breaking the iKennicuttI ( 1 1 9981) 
relation requires ^ lO" Mq kpc^^ gas (and ultimately stellar) sur- 
face densities, over spatial radii of ~ 1 — 2kpc. Although this is 
indeed comparable to the highest stellar surface densities in el- 
lipticals, observed systems only reach those surface densities at 
< lOOpc scales; not at kpc scales (where the typical surface den- 
sity in a massive elliptical is ~ IO'^Mq kpc~^). Also, given the 
characteristic SMG lifetime of ~ 10* yr, from observational con- 
straints, it is straightforward to estimate the total amount of stellar 
mass that would be formed at or above these surface densities, from 
SMGs over the redshift range z ~ 2 — 4, a nd the number comes to 
~ lO'*M0Mpc~^[E. > 10" M0 kpc~^1. iHopkins et alj l l2009j) 
calculate the actual observed stellar mass density in all galaxies 
above such a threshold, and find that it is « 1.5 x 10'' Mq Mpc~^, 
well below the range of uncertainties owing to e.g. stellar mass loss 
and the somewhat uncertain values above. 



3.3.4 The Burst Luminosity Density 

Figure [8] integrates the luminosity functions shown in Figure |7] to 
give the total IR luminosity density, and corresponding total SFR 
density, owing to bursts, as a function of redshift. Note that many 
of the uncertainties that might affect our other results are integrated 
out here (e.g. the exact lightcurve shape, or Kennicutt-Schmidt 
slope). We compare with the total SFR density from a number 
of observations in IR, rad io, and UV wavelengths as compiled in 
iHopkins & BeaconJ ( l200i) . as well as that inferred at high redshifts 
from observations of the Lyman-a forest in quasar spectra and 
implie d ionizing background, compiled in iFaucher-Giguere et al.l 
1 I2OO8I) . As expected from Figure |7] the predicted SFR density ow- 
ing to bursts is well below the total SFR density observed. This 
clearly demonstrates that bursts do not, and cannot, dominate the 
total SFR density. 

We compare with a number of observations that attempt to es- 
timate the luminosity or SFR density induced by mergers. First, we 
stress that this is different from the total SFR density in objects 
identified as mergers. A merger can take ^ 2 Gyr to complete, but 
httte. SFR over most of that time will simply reflect that of the qui- 
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Figure 8. Total IR luminosity density (and corresponding SFR density) as a function of redshift. We show the inferred contribu tion owing to bursts (solid 
black line) from the observed z = spheroid properties. We compare the observational compilation from lHookins & Beaconj j2006l) for the total IR luminosity 
density /SFR density (green diamonds), and at high redshifts the SFR density infeiTed from Lyman-o forest measurements in lFaucher-G iguere et al.' ('200^ 
green circles). We also co mpare observed estim ates of t he SFR density induced in mergers (blu e points), fr om Jogee et al. (20oJ triangles), Sanders & MirabeJ 
il99d inverted triangle), iRobaina et al.l J2OO9I circle), ISrinchmann et al.l Jl998l squares), and lBell et alj {20051 star). The burst component of star formation 
is a small fraction ~ 5 — 10% of the total IR luminosity density (most stars are formed quiescently in disks and are then violently relaxed into spheroids). 
This agrees well with observational estimates of the fraction of star formation owing specifically to mergers, as predicted if these mergers induce the angular' 
momentum loss that drives gas dissipation in forming spheroids. 



escent disks, with at most a modest enhancement. The burst itself, 
where gas dissipation drives inflows into the center of the merger 
remnant, enhancing the SFR significantly, has a duration of only 
~ 0. 1 Gyr. Nevertheless, using statistical samples of mergers and 
non-merging systems, and comparing their SFR distributions, vari- 
ous authors have attempted such a comparison.^ 

Typically, in these cases, the SFR density of some merger sam- 
ple (identified in a similar manner) is considered, but only after 
subtracting away the expected contribution from quiescent star for- 
mation. In general, this is accomplished via comparison to some 
control sample of star-form ing galaxies with sim ilar stellar masses 



and at comparable redshifts. iRobaina et al.l l 2009l) attempt this from 
a pair-selected sample at z ~ 0.4 — 0. 8; Ijogee et alj ( I2OO9I) per- 
form such an estimate from morphologically selected galaxies at 



^ Note that most of these authors actually measure the fraction of the 
SFR density in or induced by mergers, not the absolute value. We con- 
vert this to an absolute density by rescaling with the observed total SFR 
density at the same redshi ft from the best-fit observed trend presented in 
iHopkins & Beacod iioO^) . 



z ~ 0.4 — 1. We also compare with somewhat less well-defined, 
but similar samples that estimate the total amount of star forma- 
tion in observationally identified ongoing mergers or recent (mor- 
phol ogically disturbed) merge r re mnants. We compi le observations 
from lBrinchmann et alj ( Il998h and lBell et all ( l2005h . who estimate 
this quantity in morphologically-selected objects at z ~ — 1.5. We 
perform a similar exercise at low redshift using the fraction of late- 
stage major merger system s as a function of IR luminosity from 
ISanders & Mirabell \\99t ) , toge ther with the IR lumi nosity func- 



tions from 



Saunders et al 



( I1990I) and lYunetal1 ( l200lh . to estimate 
the fraction of the luminosity density in mergers. Because the selec- 
tion of ongoing mergers is somewhat strict in these samples (isolat- 
ing near-peak times), they are not very different from the estimates 
of merger-induced star formation, but formally should still be con- 
sidered upper limits. 

Our prediction agrees well with these observations, over the 
range z = — 1 . In fact, given the (probable factor ~ 2) system- 
atic uncertainties involved in both, the level of agreement is sur- 
prising, and may be somewhat coincidental. But, in general, both 
direct estimates and our indirect constraint imply that bursts consti- 
tute ~ 5 — 10% of the SFR density, and that this is not dramatically 



© 0000 RAS, MNRAS 000, 000-000 



1 6 Hopkins et al. 



larger at high redshifts. The agreement here also provides further 
evidence that most of the bright, massive bursts of interest here are 
really merger-induced, since that is specifically what is observed. If 
there were a large non-merger population, the observations should 
be much lower than our prediction. 

We have also compared both to the predictions of the mos t 
recent generation of cosmologi cal models. ISom erville et al.' (2008^ 
use semi-analytic models, and lHopkins et al. _2009g) em- 

ploy semi-empirical (halo-occupation based models), both com- 
bined with the predictions from high-resolution galaxy merger 
simulations jCox et alJboOSi : iHopkins et"al]|2009cl) to predict the 
merger-induced burst S FR and luminos i ty dens ity as a function of 
redshift. Most recently, iHopkins et al. 1 ( l2009kl) present a revised 
version of these semi-empirical models based on larger suites of 
simulations and improved halo occupation constraints and model- 
ing, and specifically present fits to the merger-induced burst pop- 
ulation. These are compared to our predicted luminosity functions 
in Figure |7] Comparing all of these results to our predictions in 
Figure [S] we find good agreement (unsurprising, given that the 
predicted luminosity functions in Figure |7] also agree well). The 
important point is that most recent models also predict that only 
~ 5 — 10% of the SFR density comes from mergers, at all redshifts. 

Finally, we hesitate to extrapolate our results beyond z ~ 2 — 3, 
for two primary reasons. First, because the Hubble time is short at 
higher redshifts, the predictions become sensitive to ~ 0. 1 — 1 Gyr 
systematic uncertainties in stellar population age estimates, which 
are known to be problematic particularly in the old-age regime. 
Second, the higher in redshift we consider, the more likely we are 
to see effects of differences between star formation and assembly 
histories, even for the bursts themselves (where, at least at low red- 
shifts, observations suggest the two are tightly coupled; i.e. bursts 
tend to be dominated by a single most massive burst, rather than 
many smaller sub-components). We therefore caution against the 
over-interpretation or extrapolation of our results at z > 4. 



4 DISCUSSION 

It has long been believed that the centers of galaxy spheroids 
must be formed in dissipational starburst events, such that gas in 
the outer parts of a galaxy must lose angular momentum rapidly, 
and fall in on roughly a dynamical time to the galaxy center. Re- 
cently, observations have shown that it is possible to robustly sep- 
arate the "burst" component of galaxy profiles from the outer, vi- 
olently relaxed component owing to the scattering of progenitor 
galaxy stars formed over earlier, more extended periods. As de- 
tailed, high-resolution observations of e.g. spheroid stellar mass 
profiles, shapes, kinematics, and stellar population properties im- 
prove, such decompositions become increasingly robust and appli- 
cable to a wide range of systems. 

In this paper, we use these observations as a novel, inde- 
pendent constraint on the nature of galactic starbursts. We stress 
that by "starburst" here, we do not mean simply any high-SFR 
system; rather, we refer specifically to star formation in central, 
usually sub-kpc scale bursts owing to large central gas concentra- 
tions driven by angular momentum loss as described above. These 
are commonly associated with galaxy-galaxy mergers, which have 
long been known to efficiently drive starbur s ts and violent relax- 
ation (e.g. Lvnden-Bell 1 9671; Toomrd 1977 ; Barnes & Hernguis^ 



Regardless of origin, we can robustly identify the starburst 
relics and stellar mass surface density profiles in well-observed 
spheroids at z = 0. Thi s methodology and tests of it s accuracy 
have been discussed in IHopkins et all ( l2008cl l2009bl H. l2008al) . 
This alone is a powerful integral constraint on the star formation 
histories of these systems. But we show that they can in fact be 
used to provide more detailed information, by allowing for the in- 
version and recovery of the star formation history of each burst. We 
can thus use local observations to recover the full time-dependent 
and scale-dependent star formation history of each burst, if we cou- 
ple these observed constraints to simple, well-motivated assump- 
tions. Namely, that some form of the Kennicutt-Schmidt relation 
(between star formation rate and gas surface density) applied in the 
formation of these stars, and that they formed in dissipational (i.e. 
rapid and initially gas-rich, on these scales) events. Both assump- 
tions are directly motivated by observations, but we also consider 
how uncertainties in them translate into uncertainties in the result- 
ing constraints. 

Performing this exercise, we recover a large number of empir- 
ically determined parameters of starbursts, and quantify how they 
scale as a function of mass and other properties. This includes the 
total mass of the starburst and the time-dependent SFR. In gen- 
eral, we find that the constrained starbursts can be reasonably well- 
approximated by a simple power-law behavior (Equation O, ris- 
ing and decaying from a characteristic maximum SFR with char- 
acteristic half-life ?burst and a typical late-time power-law slope of 
M, cc (?/fburst)~'^ where 2.0 -3.0. 

The characteristic burst mass Mbunt is ~ 10% of the to- 
tal spheroid mass, but it scales weakly with galaxy mass in 
a manner similar to how disk galaxy gas fractions scale with 
their stellar masses (expected if disks are the pre-merger pro- 
genitors of spheroids), Mbu,,,i fa 1/(1 + [M./IO^ '^M©]"*). De- 
tailed discussion of this trend, and its consequences for the 
global structu r e and kinematics of spheroids, are presented in 
IHopkins et al] bOOSal) . However, it already makes it clear that 
bursts should not dominate the SFR density. They are a small 
fraction ~ 10% of all stars in spheroids (let alone all stars). That 
does not mean that bursts are unimportant, however. It is clear 



that they control man y of the properties of g alaxies dCox et al.l 
2006l;lRobertsonetal. 2006; Naab et al. 2006; Onorbe et al.ll2006l; 



Ciotti et al.ll2007l;|je^eit et al..,2007. ,2009i ; .Covington et al.ll2008l; 
Hopkins et al. 2009bl le ). and they can account for the short-lived. 



ll99lL Il996l ; lBameslll998l ; iMihos & Hemauistlll994 Il996l) . But 

they could also owe to any other sufficiently violent process; e.g. 
gas inflows owing to disk bars or during dissipational collapse. 



highest-SFR systems in the Universe. Moreover, the scatter in burst 
mass is significant, ~ 0.3 — 0.4dex, which is critical for explain- 
ing the most extreme starbursts in the Universe, which require both 
large absolute galaxy masses and gas fractions to reach the very 
high burst masses required. 

The typical starburst timescale im plied from t he co mbination 
of observed surface densities and the iKennicuttl l ll998h law is a 
nearly galaxy-mass independent fburst ~ 10^ yr. Both the value of 
this timescale, and the weak scaling with galaxy and/or burst mass, 
agree well with the dynamical times in the central ~kpc of galaxies. 
As above, though, there is considerable scatter of order ~ 0.2 dex. 
Given such a short starburst timescale, and the relatively small to- 
tal mass fractions involved in starbursts, it naturally follows that 
starbursts will represent only a small fraction of the star-forming 
galaxies at any stellar mass, at a particular instant. Given that the 
average number of bursts per galaxy is not large, the duty cycle of 
bursts should be ~ fbursi/fHubbic, or ---^ 1 — 5% from z = — 3. Indeed, 
observations have shown that most galaxies at these redshifts lie on 
a normal star-forming seque nce, without a large 1 a scatter from 
e.g. merger-induced bursts jNoeske et al.ll2007bl lal IPapovich et al.l 
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l200dlBellet"ai]|2005h . The small duty cycle here means that even 
if the burst SFR enhancement is large, it will not violate these con- 
straints (appearing only at the ~ 2 — 3 a level in the wings of the 
SFR distribution at a given mass). 

These conclusions are supported by independent evidence 
from observational stellar population synthesis studies. Specific 
comparisons to the o bjects considered here, where a vailable, are 
prese nted in detail in iHopkins et al.l l l2009bl l3) and jposter et al.l 
I2OO9I) . It is well-established that constraints from abundances re- 
quire the central portions of spheroids be formed in a similar, 
short timescale. And detailed decomposition of stellar populations 
into burst plus older stellar populations have yielded consistent 
results for the typical bur s t fractions and sizes ( iTitus etalJI 19971 : 



ISchweizer & Seitzeilll998l : iReichardt et al.ll200ll : lMicharj|2006h . 

As should be expected from the generic behavior above, bursts 
peak at SFRs of ~ Mbursi /'burst, which follows Mbuist (and hence total 
spheroid mass) in a close-to-linear relation. The most massive local 
ellipticals - especially those with total stellar masses of > 10'^ Mq 
- had extreme peak SFRs of > lOOOMQyr"'. Thus, it is at least 
possible that some local systems reached the highest SFRs inferred 
for massive, high-redshift starburst galax ies jPapovich et"ai]|2005l : 
IChapman erauboOSi : IWalter et al.ll2009l) . More moderate, but still 
massive ellipticals with M* ~ 1 — 5 x lO" Mq, reached a range of 
peak SFRs from 30 — 500 Mq yr~ ' , corresponding to their form- 
ing fractions from ~ 5 — 20% of their total masses in starbursts. 
These match well with the observed SFRs in more typical, local 
and z ~ 1 starbursts (in ~ L, galaxies), which have been specif- 
ically associated with mergers driving gas to galaxy centers, and 
forming the appropriate nuclear mass concentrations to explain e.g. 
elliptical kine matics, sizes, phase space densities, and fundamental 
plane scalings jCqx et alj|200d iNaab et alj|2006l: [Robertson et al.l 



bOOdlJesseit et alJl2009l : lHopkins et al.ll2008jall20093r 

We similarly quantify starburst spatial sizes as a function of 
their mass, peak SFR, and time. Starburst sizes scale with starburst 
mass in a similar fashion as the spheroid mass-size relation, but 
are smaller than their host spheroids by a fraction similar to their 
mass fraction. The most massive starbursts reach half-SFR (i.e. 
half-light, in IR or mm wavelengths) size scales of ~ 1 — 5 kpc. 
For typical spatial distributions, this implies total spatial extents of 
strong emission up to ~ lOkpc. These size scales also agree well 
with obs ervations of the most massive high-redshift starbursting 
system s jYounger et alj2008l : lTacconi et alj2006l : ISchinnerer et al.l 
I2OO8I) . and of massive, compact ellipticals formed at high redshift, 
believed to be the relics of such starbursts (with, at that time, lit- 
tle envelope of dissipationless, low-density material yet accreted ; 
van Dokkum et all l2008l: Icimatti et all 120081 : iTYuiillo et alJ[20o3 : 



Bezanson et al.ll2009l : lHopkins et al.ll2009il) . 



For more typical, ~ L* starbursts, sizes range from 
0.1 — Ikpc, also similar to those observ e d in merging systems 
I Scoville et al.l Il98d ISan ders et al.l 1 19881 : Iffibbard & Yw 



I999I: 



iTacconi et al.ll2002l : lLaine et al.ll2003bl : iRothberg & Josepribooj) . 
The size difference between these and the most extreme objects 
follows from the much larger gas supply involved - there is no 
dramatic difference in the relic starburst mass distribution shapes. 
Some claims have been made that the large sizes of high-redshift 
starbursts could imply that they are not scaled up analogues of lo- 
cal extreme starbursts; but we find here that they correspond nat- 
urally. Scaling up a starbursting system in the starburst mass frac- 
tion will not preserve spatial size, but rather will scale along the 
starburst size-mass relation here, which appears to be smooth and 
continuous from the smallest starbursts with masses lO' Mq to 
the largest with masses > 10"Mq. The origin of the size-mass 



scaling is of considerable physical interest. It has been proposed 
that in such starbursts the Eddington limit from radiation pressure 
on dusty gas sets a universal maximum c entral surface density , over 
all mass scales, from which this follows ( iHopkins et alj2009j) . The 
important constraint, from our analysis, is that there is no discon- 
tinuity, and we provide the scaling relations that any such model 
must satisfy. 

Combining these constraints with observational measure- 
ments of the nuclear stellar population ages of these systems - 
i.e. the distribution of times when these bursts occurred - we show 
that it is possible to re-construct the dissipational burst contribution 
to the distribution of SFRs and IR luminosity functions and lumi- 
nosity density of the Universe. We show that the burst luminosity 
functions agree well with the observed IR LPs at the brightest lu- 
minosities, at redshifts z ~ — 2. At low luminosities, however, 
bursts are always unimportant, as expected from their short duty 
cycles, noted above. Although the burst luminosity functions rise 
with redshift, they always represent low space densities, and the 
overall LF evolves rapidly. As such, the transition luminosity above 
which bursts dominate the IR LFs and SFR distributions increases 
with redshift from the ULIRG threshold at z ~ to Hyper-LIRG 
thresholds at z ~ 2. This appears to agree well with recent esti- 
mates of the transition between normal star formation and merg- 
ers, along the observed lumino sity functions. Systematic morpho- 
logical studies at low redshifts jSanders & Mirabelll 19961) yield the 
conventional wisdom that - locally - the brightest LIRGs and es- 
sentially all ULIRGs are merging systems (see also references in 
§ [T). At high redshifts, sim ilar studies have now been performed 
(see e.g. lTacconi et alj2008l and references therein). They too find 
that the brightest sources are almost exclusively mergers, but with a 
transition point an order-of-magnitude higher in luminosity. Other 
morphological studies at int ermediate redshifts z ^ 0.4 — 1.4 have 
reached similar conclusions ( [Bridge et al.ll2007l) . 

Integrating these luminosity functions, we find the burst con- 
tribution to the SFR and IR luminosity densities of the Universe, 
and show that it is small at all redshifts, rising from ~ 1 — 5% at 
z ~ to a roughly constant ~ 4 — 10% at z > 1. This agrees well 
with recent attempts to estimate the contribution to the SFR den- 
sity at z = — 1 specifically from merger-induce d starbursts, using 
either pair or morph ologically-selected samples jjogee et^l2009l ; 
iRobaina et al.ll20()9l) .^ Given the completely independent nature of 
the constraints, and significant uncertainties involved in both, the 
agreement is good. The small value is what is expected, given our 
previous determination that the typical burst mass is just ~ 10% 
in ~ L* spheroids (the galaxies that dominate the stellar mass den- 
sity). But it clearly rules out merger-induced bursts driving the SFR 
density evolution of the Universe. 

At the highest redshifts z > 2, we can put strict upper lim- 
its on starburst intensities, based on the maximum stellar mass re- 
maining at high densities at z = 0, and find some tension between 
these and estimated nu mber counts of sub-millimeter galaxies from 
IChapman et al.l l l2005h . This implies that some change may be nec- 
essary in either the number counts themselves, the bolometric cor- 
rections used to convert these observations to total IR luminosities, 



* Note that it is important here to distinguish estimates of the SFR induced 
by mergers, i.e. that above what some control population would exhibit, 
from that simply in ongoing/identifiable mergers (since many criteria iden- 
tify mergers for a timescale ~Gyr, much longer than the burst timescale). 
For example, a merger fraction of 10% would imply at least 10% of star 
formation in ongoing mergers, even if there were no starbursts and those 
systems were only forming stars at the same rate as they would in isolation. 
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or the stellar IMF used to convert between SFR and IR luminosity. 
However, the observations remain considerably uncertain, with the 
bolometric corrections relying sensitively on assumed dust temper- 
atures, and th e number counts subjec t to significant cosmic vari- 
ance (see e.g. lAustermann et alJIlOOgh . More observations, at new 
wavelengths and in larger, independent fields, are needed to resolve 
these discrepancies. 

We compare our constraints on these histories with recent pre- 
dictions from galaxy formation models and simulations, and find 
reasonably good agreement. Both exhibit similar tension with es- 
timates of the high-redshift, extremely luminous number densities. 
However, the models are able to match the inferred number den- 
sities presented here without any change to the stellar IMF or re- 
quiring other exotic physics. The systematic uncertainties in the 
models, especially at high redshift, are large, however, so the em- 
pirical constraints presented here provide a powerful new means of 
constraining the models and their input parameters. 

For example, if models are constrained via e.g. a halo occupa- 
tion approach or otherwise, so as to match the observed merger 
fractions of galaxies as a function of redshift, then these num- 
bers become relati v ely large (> 10%) at r edshifts z > 2 (see e.g . 
Bundv et alj l2()09l: IConselice et alj [200^; iKartaltepe et alj [2007 



Lin et alj200l:lBluck et alj2008l ; lHoDkins et al.l2009ll : |jogee et al.l 



2009l : lBridgeetalJl2009l) . The most common assumption m many 
analytic and semi-analytic models is that, in such a major merger, 
the entire galaxy gas supply is channeled into the starburst. This, 
however, coupled with the high observed merger fractions (and im- 
plied merger rates), would lead to an over-prediction of the burst 
contribution to the SFR density, relative to our const raints here. For 
exampl e, the predictions of such a model are given in lHopkins et al.l 
( l2006hl) . and rise t o ~ 20 - 50% of the S FR density at z > 1 - 2. 
Recently, however, iHopkins et al] j2009d) have pointed out that the 
processes that lead to angular momentum loss by the gas - driving 
bursts in the first place - rely on the stellar component of merging 
systems, and so become less efficient as gas fractions increase. In- 
cluding these physics in the models (as the current generation to 
which we compare does) leads to an asymptotic maximum con- 
tribution to the SFR density similar to that found here, while still 
giving a good match to observed merger fractions. 

Despite the simple nature of the assumptions involved in our 
analysis, we show in high-resolution hydrodynamic simulations 
that they work well in allowing us to recover the star formation 
histories of bursts (even where the simulated system is much more 
complex). Especially in a statistical sense, our numerical experi- 
ments suggest that this approach is robust to a variety of detailed 
deviations from our idealized assumptions. Testing this (even in 
a few objects) via high-resolution reconstruction of stellar pop- 
ulations in bursts, directly from their observed spatially-resolved 
SEDs, would provide a powerful test of this. Insofar as we have 
compared with different simulations and models, the dominant un- 
certainty in our analysis is the nature - in particular the slope, rel- 
evant for t he ext rapolation to high gas surface densities - of the 
iKennicut 3l ll998l) law. We have considered a range of slopes sug- 
gested by different observations. Although they yield similar be- 
havior at moderate and low SFRs, after approximately ~ 10* yr 
from the beginning of a burst, the predicted behavior at very early 
times in the bursts is quite different. A larger [Kennicutti ( [l99S ) re- 
lation index implies higher peak SFRs and more sharply peaked 
starbursts, increasing the predicted number counts of the most lu- 
minous sources and making it relatively more easy for models to 
account for the most extreme star-forming systems. It therefore re- 
mains of considerable importance for observations to probe both 



gas density measurements and full SFR indicators in extreme sys- 
tems, at low and high redshifts. 

Finally, we note that our analysis is only possible because, 
as indicated by basic dynamics, simulations, and observations of 
stellar populations, starburst components of spheroids were formed 
in dissipational (i.e. gas-rich, at their centers), rapid star forming 
events. The dissipationless "envelopes" surrounding the central, 
dense components in spheroids were not formed in such a manner 
(again indicated by both their structural and kinematic properties, 
simulations of their formation, and stellar population observations). 
Rather, they represent the debris of stars from disks which were 
formed pre-merger, and assembled (and violently relaxed) dissipa- 
tionlessly. These stars were formed over extended periods of time, 
with new gas accr etion onto the disk fuelin g new or continuous 
star formation (e.g. iKeres & Hernauislll2009l) , as opposed to a sin- 
gle massive inflow. Especially in the most massive systems, they 
are also assembled from multiple systems, via e.g. minor mergers 
contributing tidal material to the extended "wings" well-known in 
massive galaxies. As such, there is not a straightforward means to 
invert their surface stellar mass density profiles to obtain their star 
formation history. Such constraints will depend on other methods, 
such as e.g. direct stellar population analysis. It should be borne in 
mind that this represents ~ 90% of the mass in most spheroids, and 
so understanding star formation in disks remains critical to under- 
standing the origin of stellar populations in spheroids. 
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Table 3. Fits to IR LF of Bursts 

LoCzbALo)" L'" L"'- 00" <t>"' 



Double Power-Law Fit - Standard Model "' 

11.87(0.09) 5.08(0.43) -7.29(0.48) -6.29(0.17) 3.76(0.85) -2.48(0.94) 0.98(0.09) 0.39(0.49) -3.33(0.62) 2.68(0.16) 7.13(0.82) -11.38(0.94) 



Double Power-Law Fit - Steeper Kennicutt-Schmidt Index " 
11.89(0.09) 6.75(0.44) -9.08(0.51) -6.27(0.15) 1.47(0.78) -0.15(0.90) 0.87(0.07) 0.63(0.40) -2.63(0.49) 2.29(0.12) 8.11(0.70) -12.32(0.84) 



Parameters are given for best fit to the redshift-dependent form of the IR LF (Equations [16117): ^{z) = <j), /[{L/U)" + {L/Uyi]. 

° '' Break luminosity L» as a function of redshift, per Equation[T7] log {L, /Lq} = Lo + L' ^ + L"5~, where ^ = log (1 +z) 
''"^ Normalization 0* : log {0* /Mpc"' log" ' Lir} = 0o + 0' ^ + 0" 
s ' Faint-end slope o: a = oq + c*' 5 + c*" 
J"' Bright-end slope 7: 7 = 7o + 7' C + 7" 

'" Luminosity functions assuming the low-redshift fitted slope for the Kennicutt-Schmidt relation between star formation and gas surface density, = 1 .4 
" Same, assuming instead the steeper slope suggested by high-redshift observations, hk = 1.7 

All fits are based on z = — 4 results. Extrapolation beyond these redshifts should be performed with caution. Statistical uncertainties in the fitted parameters ai'e given in 
parentheses. 



